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Abstract 

We study how much gluon shadowing can be perturbatively generated 
through the modified QCD evolution in heavy nuclei. The evolution of 
small-x gluons is investigated within the semiclassical approximation. The 
method of characteristics is used to evaluate the shadowed distributions 
in low-Q and small- a; region. In solving the modified evolution equation, 
we model in simultaneously fusions from independent constituents and 
from the same constituent, both in a proton and in a large loosely bound 
nucleus of A ~ 200. In addition to the actual distributions at small x, we 
study the ratios of the distributions at an initial scale Qo = 2 GeV, and 
show that a strong nuclear shadowing can follow from the modified QCD 
evolution. 
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The semihard gluonic subprocesses are expected to play an essential role in the 
formation of high energy densities in heavy ion collisions at collider energies |1|, 0, ^, ^. 
However, there are many theoretical uncertainties in modeling these QCD processes. 
One of the major ones, nuclear gluon shadowing, comes from the unknown initial 
gluon distributions at small x. Unlike for the quark and antiquark distributions, 
there are no direct experimental data for the gluons in nuclei. Getting theoretical 
control over the nuclear gluon shadowing is therefore a very urgent and important 
issue. The purpose of this Letter is to study how much gluonic shadowing is generated 
perturbatively through the modified QCD evolution [§, ^ in heavy nuclei. 

"Shadowing" in the context of the deep inelastic ZA-scattering refers to the mea- 
sured depletion of the nuclear structure function at small xbj, as compared to F2 
of unbound nucleons 0. The same kind of depletion at small x is expected to happen 
also in the nuclear gluon distributions. During the recent years there have been many 
efforts to explain the measured nuclear shadowing of quarks and antiquarks |]8[|-|[T5[| 
but for gluons the situation is still inconclusive. Once the nuclear parton distributions 
are known at an initial scale Qq, the QCD-evolution to larger Q can be computed 
[0, O, O. The problem is how to get input, theoretically or experimentally, for the 



nuclear gluon distributions at Qq, and to understand the reliability of QCD-evolution 
for the proper range of x- and Q- values. 

Shadowing-phenomenon is also predicted to happen in protons. In this case, 
"shadowing" refers to the depletion of the actual parton distributions, and is caused 
by the fusions of overcrowding gluons at very small x. This mechanism proceeds 
through perturbative QCD-evolution as formulated in |^, ^. It has been shown by 
Collins and Kwieciiiski that the singular gluon distributions actually saturate due to 



the fusions [|^. 

In this Letter our basic idea is quite straightforward. We first compute the gluon 
shadowing in a proton, by using techniques introduced in |jl8[ for solving the small- 
X evolution including gluon recombination. Then we apply the same mechanism of 
recombining gluons to heavy nuclei and study to what extent nuclear shadowing is 
generated perturbatively through the QCD-evolution at Qo = 2 GeV. 

At small values of x, leading order QCD evolution equation predicts that the num- 
ber of gluons becomes extremely large. It has been known p, Q that for sufficiently 
small values of x and/or of Q^, the total transverse area occupied by the gluons will 
be larger than the transverse area of a hadron, so that the interaction between gluons 
can no longer be neglected. Such gluon recombination results in a modification of 
the QCD evolution equations. In the limit of small-x the modified QCD evolution 
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equation can be cast in the form p, ^ 

dydtGiy, t) = cG{y, t) - X exp(-t - e*) [G{y, t)]^ (1) 

where y = ln(l/x), t = ln[ln(QVA^cD)]> Giv, t) = ^di^, Q^) and c = 12/{ll-2Nf/3) 
with Nf the number of quark flavors. 

Strength of the gluon recombination is controlled by the factor A, originating from 
two possible sources. The two fusing gluon ladders, which couple 4 gluons to 2 gluons, 
can arise either from independent constituents of proton/nucleus or from the same 
one, as discussed in |]^, |TB|, |2T|. We will refer to the former "independent" and 

to the latter as "non-independent" fusion. Since recombinations from both sources 
happen simultaneously, we divide the parameter A into two parts: 

A = Ai + An, (2) 

where Ai corresponds to the independent recombination and An to the non-independent 
one. 

Let us first study the two sources of recombination within the models for two- 
gluon densities given in ref. In a proton, the strength of the independent fusion 
then takes the form 

2 1 TT^C^ 

= 3^ ■ 2A|^' 
where i?p ~ 1 fm is the radius of a proton. 

The magnitude of the non-independent fusion of the gluon ladders can be esti- 
mated as 

16 1 TT^c' 

"^8l7r(2/Q,)2'2A2jCD' ^ ^ 

where we have made a simplification by fixing the initial x of the valence quark to 
~ 1. We also approximate the scale of the initial valence quark by ~ 2 GeV. 
Let us then consider a large loosely bound nucleus. Naturally, both types of 
fusions are still there but only for the independent one an A^/^-scaling arises. In this 
case 

^ 9 A TT^C^ 

= 8^ ■ 2A|^' 

where the nucleus is taken to be a sphere with a sharp surface at Ra = 1.12y4^/^ fm. 
The strength of the non-independent fusion remains the same as in the case of a free 
proton: Aj^ = An. 

It is interesting to notice how the relative contributions of the two types of recom- 
bination will change when going from a proton to a nucleus of A ~ 200: Xu/Xi ~ 7.6 
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and XfJXf ^ 1.0. Thus the non-independent fusion is clearly dominant in a free 
proton while in a large nucleus the contributions from both types are of the same 
order. As a result, parton recombination is strongly enhanced in a heavy nucleus. 

In order to solve Eq. |l| exactly by integration, one would need the initial distri- 
bution either at fixed or to o.'^^d the derivatives along a boundary line (y, to) or 
{yo,t), respectively. However, since the expression for the non-linear term in Eq. |l]is 
not valid for the regions where x is large, or where both x and Q are very small, the 
natural boundary condition at a; = 1 (or yo = 0) is not suitable here. In addition, 
since we do not have sufficient information on other boundary lines, we cannot solve 
Eq. by direct integration. Instead, with the semiclassical approximation , we are 



going to adopt the idea introduced in [18| to use the method of characteristics, so 



that we can avoid the region III (see discussion later). 

The semiclassical approximation corresponds to neglecting the second order deriva- 
tive term, dydt ln{G), which leaves us with the evolution equation as 

dyz{y,t)dtz{y,t) = c - Xexp[-t - + z{y,t)], (6) 

where z{y,t) = \n[G{y,t)]. The above equation can then be cast and solved in the 
form of a set of characteristic equations as shown in detail in WE\. 



The evolution of gluon distribution in a proton and in a nucleus is similar, so let us 
first consider the general idea, as illustrated in Fig. 1. We divide the (y, t)-plane into 
three regions. In region I for y < yo{x > xq), we expect the traditional, non-corrected 
Altarelli-Parisi (AP) evolution [|l9l to hold down to scales Q ~ 1 GeV. In region II 
with y > yo{x < xq), the evolution of the gluon distribution is then approximately 
given by Eq. The initial values of the gluon distribution and its t- and ^/-derivatives 
are determined numerically at (?/o,tmin < < ^o) from the singular gluon distribution 
of the CTEQ-collaboration [Q, which is obtained with a lower initial Q- value to 
provide the necessary t-dependence of the boundary condition at yo, and a larger 
x-cut to ensure the validity of AP evolution in region /. The evolution in the region 
// is finally terminated at t = to, corresponding to Qo = 2 GeV. Notice that the 
characteristics approach the to-line from below; the lowest scale we have to go down 
to is about Q = 1.25 GeV, corresponding to tmin- At these scales QCD perturbation 
theory should be still valid in region /. In region /// with extremely large y (small 
x) and/or small t we do not expect our analysis to be valid anymore, since the higher 
order terms in the evolution equation will become important. 

Before performing the actual evolutions, we have to consider how to choose the 
boundary yo{= — In(xo)) for a proton and a nucleus, and how to conserve momentum. 
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For a proton, we assume that the recombinations start to be effective at a; ~ xq 



0.01, which is consistent with ||T8| , We can use the resuhs from global fitting, like 
CTEQ, to constrain xq. In fact, we will see that with xq = 0.01 the shadowed gluons 
deviate considerably from the CTEQ gluons only after x < 0.001, so the choice for 
Xq seems to be reasonable, and we do not expect the results to be very sensitive to 
small changes of xq- 

As explained above, the gluon recombination is strongly enhanced in heavy nuclei 
and it starts at somewhat larger values of x than in protons. The corresponding 
boundary line Xq for a nucleus is approximately determined by the relative magni- 
tude of the evolution terms in Eq. ^ Ga{xq) ~ G{xo)\a/ ^, so that the relative 
contribution from the gluon fusion in a nucleus is about the same as in a nucleon. 
This gives Xq ~ 0.05-0.1. This range of Xq is also supported by other studies [|1^ . 

Let the original fraction of momentum in gluons be /o = Jq dxxgcTEQix,Ql). 
In the case of a proton, shadowing in the region II changes the gluonic momentum 
typically by less than a per cent, which we can clearly neglect as a small overall 
change. 

Perturbative shadowing reduces the gluonic momentum more in a nucleus than in 
a proton. Assumed that the momentum fraction of gluons is conserved, there must 
be a corresponding enhancement in the region II. In addition to this, we also take 
into account a possible momentum transfer from quarks and antiquarks to the gluons. 
Here we consider nuclei with A ~ 200, for which we expect an overall increase in the 
fraction of the momentum, e^, to be only about 4% |P, |10|, |T^. We combine these 



two sources of the momentum flow, which results in solving iteratively from 

/ dxxg{x,Ql) +aA dxxgcTEQix,Ql) = fo{l + eA), (7) 

Jo C Jx^ 

with the condition C : gA{xQ,Ql) = aAfi'CTEQ(a^o , Qo) '^^ boundary. Typically, 
aA ~ 10% for A ~ 200. 

Let us now turn to the results, presented in Figs. 2. In Fig. 2a, nucleon and 
effective nuclear gluon distributions for a nucleus of A = 200 are compared with the 
input CTEQ gluon distribution at Qo = 2 GeV. Notice the ~20 % uncertainty in the 
nuclear case resulting from varying Xq from 0.05 to 0.1. To demonstrate the formation 
of strong perturbative nuclear shadowing, corresponding to the relative depletion of 
gluon distributions in a nucleus, we plot the ratio Ga{x,QI)/G{x,QI) in Fig. 2b. 
Notice also that as x decreases, gluon distribution in a proton increases much faster, 
or shows the sign of saturation at a much smaller x than that in a nucleus. Therefore, 
as shown in Fig. 2b the ratio saturates only when the gluons in a proton do so. Thus, 
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saturation of the perturbative nuclear shadowing reflects actually the behavior of the 
gluons in a proton. As a main result, we conclude that, due to the enhanced gluon 
recombination in a heavy nucleus, a ~50% nuclear shadowing in small-x region is 
generated perturbatively through the modified QCD evolution, accompanied by a 
~10 % antishadowing from the momentum conservation. 

As seen more clearly in Fig. 2b, with Xq = 0.1 there is a slight deviation in the 
initial derivatives of the gluon distributions determined from the region / as compared 
to what can be determined from Eq. ^ from the region //. This in turn is a reflection 
of an apparent fact that eventually one cannot apply the small-x approximation at 
too large x. Letting Xq ~ 0.1, we are really pushing the small-x evolution equation 
to its limit; surely beyond this point the Eqs. |l] and ^ cannot be applied without 
additional correction terms. However, taking the initial conditions from the "known" 
region /, as we do, should improve the analysis and reduce the uncertainty in the 
small-x region. From Fig. 2b it is seen that we cannot make conclusive claims about 
the "beginning" of nuclear shadowing. However, since the result with Xq = 0.1 does 
not differ considerably from the result with Xq = 0.05, we believe our result shows 
the correct order of magnitude of the perturbative shadowing at very small x. 

It is clear that the absolute strengths of the A's depend on the models assumed 
for the two-gluon densities in a proton and in a nucleus. However, we do not expect 
our qualitative results for the nuclear shadowing to change very much with different 
details. One may also question what happens to nuclear shadowing, if the initial 
gluons diverge more strongly(weakly) when a; than CTEQ gluon distribution 
used here. In that case, the recombinations would be enhanced(suppressed) both 
in a proton and in a heavy nucleus. However, the result for nuclear shadowing, 
Ga{x,QI)/G{x,QI), is not extremely sensitive to the small-x behavior of the input 
gluon distribution because we use only the part with x > 0.01, which has been 
relatively well-tested experimentally. This question will be studied in more detail 
elsewhere . 



We would like to comment briefly on the general consequences of our result for 
perturbative nuclear shadowing. The semihard processes with typical scales Q ~ 
a few GeV involve x ~ Qj^. In heavy ion collisions at ^/s = 200 GeV, the x's 
in the semihard processes will be larger than 0.01, so these processes rather probe 
the onset of perturbative nuclear shadowing than the region of saturation. On the 
contrary, in collisions with ^/s in the TeV range, the semihard processes will happen 
at x's typically smaller than 10"'^, and are therefore affected considerably more by 
the perturbative shadowing. Examples of the possible effects on minijet production 
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can be found in 0, |17|, Other processes clearly suppressed by the nuclear gluon 
shadowing at very high energies are heavy quark and their bound state production. 
Also the production of total transverse energy and energy density will be suppressed, 
as compared to the predictions with non-shadowed gluons p|, Through this, the 
thermalization of the possibly formed quark-gluon plasma is also slowed down, in 
which case the thermal electromagnetic signals are suppressed. In order to make 
precise predictions for these processes, nuclear shadowing has to be studied at scales 
Q > 2 GeV. The scale dependence of nuclear gluon shadowing is an interesting 
question to which we will return in the future [^ . 

To conclude, we have considered the perturbative aspects of the nuclear modifica- 
tions to the gluon distributions. As we have shown here, a strong nuclear shadowing 
is generated through the modified QCD evolution, and it may well be the dominant 
mechanism for the small-x modifications. We emphasize that the use of the method of 
characteristics is necessary to avoid the region III where even the modified evolution 
equation is not expected to be valid. We feel we now have more quantitative con- 
trol over the nuclear gluon distribution at small a;, based on perturbative QCD. We 
believe this study could serve as an interesting starting point for more detailed cal- 
culations of nuclear gluon shadowing and its consequences in ultra-relativistic heavy 
ion collisions. 
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Figure Captions 

Fig. 1. The evolution plane. In the region I the traditional AP-equations are 
expected to be valid. Both x- and Q-dependence in this region form the initial 
conditions for the evolution in the small-x (large y) region II . Examples of the 
characteristics of the Eq. ^ in the region // are shown. 

Fig. 2. a. The gluon distributions xq^x^Qq) at Qo = 2 GeV vs. x. The result for 
proton is labeled by xq, and the results for A ~ 200 by x^, respectively. The CTEQ 
gluon distribution [^ is labeled by "CTEQ". b. The ratio xqa^x^Qq) /xg{x,Ql) of 
the shadowed gluon distributions vs. x, demonstrating a strong perturbative nuclear 
shadowing in heavy nuclei. 
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